Pump-controlled directional light emission from random lasers 
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The angular emission pattern of a random laser is typically very irregular and difficult to tune. 
Here we show by detailed numerical calculations that one can overcome the lack of control over this 
emission pattern by active control of the spatial pump distribution. We demonstrate, in particular, 
how to obtain customized pump profiles to achieve highly directional emission. Going beyond the 
regime of strongly scattering media where localized modes with a given directionality can simply 
be selected by the pump, we present an optimization-based approach which shapes extended lasing 
modes in the weakly scattering regime according to any pre-determined emission pattern. 

PACS numbers: 42.55.Zz, 42.25.Fx 



The scattering of waves in a disordered medium is 
a ubiquitous physical process that occurs in a variety 
of different scenarios ranging from the micro- to the 
macroscale [1, 2]. Although much attention has been 
devoted to this topic, new intriguing phenomena con- 
tinue to emerge in this area of physics [3-6] . Particularly 
exciting progress towards practical applications has re- 
cently been reported from the field of disordered photon- 
ics [7]: randomly scattering materials with optical gain 
were shown to operate as lasers despite the absence of 
a cavity [8-10]. In contrast to conventional lasers, how- 
ever, such random lasers emit light at multiple frequen- 
cies and with a broad spatial emission profile [8-12]. On 
the other hand, unprecedented control of light fields in 
complex media has been achieved through optical wave 
front shaping based on spatial light modulators [13]. This 
progress has, e.g., enabled the spatial and temporal fo- 
cusing of light or the projection and retrieval of images 
through strongly scattering media [14-18]. 

In the present letter we strive to bridge the advances 
in the above two fields in order to tackle the longstanding 
goal [19-33] of making random lasers externally tunable. 
Employing a quasi-one-dimensional (ID) laser, it was re- 
cently shown numerically [31] as well as experimentally 
[32] that, using wave front shaping techniques to iter- 
atively optimize the pump profile applied to a random 
laser, single-mode lasing operation can be achieved at a 
predetermined frequency. With this approach one of the 
major deficiencies of a random laser, i.e., the lack of con- 
trol over its emission frequency could be remedied. This 
leaves random lasers with the major drawback of emit- 
ting into arbitrary directions - a property which severly 
limits their applicability in many practical circumstances. 
Ideally, of course, the control over the angular emission 
pattern could also be exerted through the applied pump 
profile which is the most easily tunable control knob in a 
random laser. Demonstrating this explicitly seems, how- 
ever, a formidable theoretical challenge. Consider, for 



^^^^^^^^ Pump 

Dielectric scatterers 




FIG. 1. (Color online) Sketch of the envisioned setup for 
realizing pump-controlled emission from a random laser: Be- 
fore being directed to the random medium, an applied pump 
beam is sent through a spatial light modulator (SLM) which 
shapes the beam following an optimization procedure. The 
targeted emission pattern can be chosen at will. The geometry 
of the disk-shaped random laser (radius R = 1 \im) consists 
of square shaped scatterers with index n s (in red), embedded 
in a uniform background medium with index nb = 1 (gray). 



comparison, the complexity involved in the demonstrated 
pump optimization for random lasers [30-33] or in the 
well-studied problem of engineering a microcavity laser 
with unidirectional emission [34-40] . In both of these two 
cases, the functions which were optimized in one or sev- 
eral sequential steps were strictly one-dimensional: For 
the random lasers this function is related to the spatial 
pump profile and for the microcavity lasers it is the cavity 
boundary. In contrast, for fully controlling the emission 
pattern of a random laser one has to deal with an es- 
sentially two-dimensional pump profile with a very large 
number of adjustable degrees of freedom. 
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We address this challenge numerically, employing stan- 
dard semiclassical laser theory at lasing threshold [41]. 
Our dielectric random laser medium consists of a pas- 
sive and active component, e = e c + Xg-> with the passive 
dielectric function e c given by a static and randomly gen- 
erated as well as immobile distribution of square-shaped 
sub- wavelength scatterers (see Fig. 1). These scatterers, 
which all have the same index of refraction n s , are sus- 
pended in a disk-shaped background material of index 
rib — 1- Both the scatterers and the background medium 
are assumed to be optically active, with a quantum sus- 
ceptibility Xg that depends on the wavenumber k [42], 

Xg = D F( X )j(k), with 7 (fc) = - 7± . . (1) 

This active component is subject to a position-dependent 
external pump source F(x) as realized, e.g., with an op- 
tical pump beam controlled by a spatial light modula- 
tor (see schematic setup shown in Fig. 1). Here, Do 
defines an overall pump strength, k a is the frequency 
of the atomic transition (we consider the speed of light 
c = 1) and 2jj_ is the width of the Lorentzian-shaped 
gain curve \j(k)\ 2 . Note, that this approach is identi- 
cal to the non-interacting form of the recently developed 
steady-state ab-initio laser theory (SALT) [12, 43, 44], 
which provides the stationary solutions of the semiclas- 
sical Maxwell-Bloch equations. We can therefore make 
use of several parts of the efficient methods developed 
for SALT in the optimization problem at hand. In par- 
ticular, we will employ the so-called threshold constant 
flux (TCF) states u n defined by the non-Hermitian eigen- 
value problem, 

{V 2 + k 2 [e c (x) + 77 n (fc)F(x)]} Mn(x, k) = . (2) 

These states n n (x, k) and the corresponding eigenvalues 
r) n (k), which we parametrize over the wavenumber k out- 
side of the disk, are self-orthogonal with respect to the 
pump profile F(x). They furthermore have the impor- 
tant property that each threshold laser mode (TLM) \D^, 
which we are interested in, is identical to one of these 
states, %(x). (Note that, throughout the paper, we will 
use overbars to denote quantities at the lasing thresh- 
old.) The corresponding required pump strength Dq and 
frequency k{ at threshold are obtained from the thresh- 
old condition 7] n (ki) = 7(A^)Z)q [44, 45]. We numeri- 
cally calculate the TLMs using a high-order finite ele- 
ment method which efficiently copes with the finite-sized 
scatterers and the spatially varying pump profile F(x) 
[46]. 

The task that we address within the above framework 
is to shape the emission pattern of the TLM which has the 
lowest (first) pump threshold Dq. Considering the scat- 
terers of our disk-shaped random medium £ c (x) as well 
as the gain curve j(k) to be immobile (as in most experi- 
ments), we only tune the applied pump profile F(x) such 



that the laser emission at threshold becomes unidirec- 
tional. To quantify the unidirectionality of a laser mode 
we evaluate its far- field profile FFP(^) for a certain pump 
configuration, which is characterized by a set of complex 
expansion coefficients {Pm,n} in a Bessel- function basis 
that respects the circular symmetry of the disk boundary 
[47]. Next, we introduce a weight V which measures the 
overlap of the far-field profile FFP(^) of this mode with 
a desired profile G(ip), 

where the denominator ensures that the weight is nor- 
malized such that V £ [0,1]. Choosing for the target 
function G(<p) a narrow Gaussian profile with a peak cen- 
tered around ip — 180° (see red curve in Fig. 2c), we can 
formulate the goal of unidirectional emission in this di- 
rection as an optimization problem for the coefficients 
{ftm,n} to be tuned such that V —> 1. 

Before going into the details of the optimization pro- 
cedure, let us first generally assess the applicability of 
our approach for a dielectric featuring the two compet- 
ing terms, e = e c + Xgi °f the passive and the active 
medium, respectively. Depending on which one of the 
two terms dominates, the structure of the lasing modes 
will be determined either by the passive medium or by 
the pump function, respectively. Correspondingly, our 
approach based on controlling the pump profile will have 
to be customized for each of these two limits. 

In the regime of strong disorder, where e c dominates 
over Xgi the disorder itself will strongly trap modes with, 
correspondingly, a high Q-factor and a low lasing thresh- 
old. The structure of these lasing modes is thus mostly 
determined by the distribution of scatterers and not by 
the pump profile. Rather than optimizing the latter, we 
can here just choose a pump profile that simply selects 
the mode with the smallest deviation from the desired 
unidirectional far-field pattern to be the first lasing mode 
[33, 48, 49]. 

The situation is very different in the opposite limit of 
weak disorder, where Xg dominates over e c . As no high- 
Q modes exist from the outset, also mode- selection does 
not work. Rather, the pump profile has to be optimized 
such as to shape the modes in a way to be explored be- 
low. Since the case of strong disorder yields directional 
emission in a straightforward way, it provides a natural 
benchmark to gauge the quality of the emission patterns 
obtained for weakly disordered samples. We will thus de- 
scribe both the weak and the strong disorder case, which 
we realize by a different contrast between the refractive 
index of the rectangular scatterers n s in our sample and 
the background index, rib = 1. 

Let us first examine the limit of strong disorder, for 
which we choose n s = 3. In this limit, the TCF 
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FIG. 2. (Color online) First threshold laser mode ^1 at 
ki = 29.82 u-m -1 of a random laser in the strongly scat- 
tering regime (see Fig. 1 with n s — 3.0). Shown are (a) 
the mode profile |^i(x)| inside the random laser, (b) the 
Gaussian-shaped pump profile that selects the mode (darker 
means more pump) and (c) the far-field profile (blue) as well 
as the desired emission profile used as a reference (red). The 
inset in (c) depicts the near field directly outside the random 
laser disk and is complementary to (a). Note that the best 
directionality of a localized mode (V = 48.19%) matches the 
desired profile only approximately. The employed gain curve 
has a center frequency at k a = 30.0 urn -1 and a width of 
7_i_ = 0.2 \irn~ 1 . For better visibility the linear color func- 
tions used in the plots differ in absolute scale. 

states linfx, k) resulting from Eq. (2) with uniform pump- 
ing, F(x) = const., are strongly localized in space [49]. 
Among these high-Q modes we also find many modes 
near the boundary of the disorder region with, corre- 
spondingly, a very narrow emission profile. To achieve 
the goal of directional emission, e.g., in the direction of 
tp = 180°, we just need to focus the pump beam on the 
mode that emits most pronouncedly in this direction to 
make it lase first [49]. We demonstrate this explicitly 
with the example shown in Fig. 2, where a directional las- 
ing mode is considered, which is the second TLM for the 
uniformly pumped device. Focusing a Gaussian shaped 
pump beam on the center of the region where this mode 
is localized, we can rearrange the modes such that this 
desired mode has the lowest lasing threshold of all. We 
have explicitly tested this pump selection procedure for 
a number of modes with different orientation of direc- 
tional emission and found it to work very well, provided 
that enough directional modes are available and the gain 
curve is broad enough to be able to select them. 

In the weak disorder limit, for which we choose n s = 
1.2, the lasing modes are extended, overlapping and their 
shape strongly depends on the pump profile F(x) [12]. 
Simple mode-selection is thus ruled out and a more so- 
phisticated optimization is required instead. We start 
with an adjustable pump profile characterized by 153 
complex expansion coefficients {Pm,n} [47]. An optimiza- 
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FIG. 3. (Color online) First threshold laser mode ^i at 
ki — 30.00 u-m -1 of a random laser in the weakly scattering 
regime (see Fig. 1 with n s — 1.2). Shown are (a) the opti- 
mized mode profile |^i(x)| inside the random laser, (b) the 
corresponding pump profile resulting from the optimization 
algorithm and (c) the far-field profile featuring a directional- 
ity of V — 96.08%. The inset in (c) shows the near field of 
the laser mode, just outside the pumped disk. Note the very 
good agreement which we find in the main panel of (c) be- 
tween the optimized emission pattern (blue) and the targeted 
Gaussian emission profile G((p) (red) centered at (p = 180° 
and with a width of 8.3°. In (d) the directionality is shown 
with respect to the number of iterations used in the optimiza- 
tion procedure. The parameters of the employed gain curve 
are k a — 30.0 u-m -1 and 7^ = 0.01 um -1 . 



tion procedure is now employed to determine these pa- 
rameters such that the far- field emission profile FFP(^) 
of the first TLM follows the desired narrow shape G(cp) 
or, equivalently, V — » 1 in Eq. (3). Note also that 
it is clearly unfeasible to find the unique global maxi- 
mum for V(^i) in the extremely large parameter space 
of the above complex coefficients. Instead, we start with 
random initial guesses for the pump profile coefficients 
{Pm,n} which we further refine locally by a gradient- 
based optimization [50]. For the local optimization pro- 
cedure one needs to evaluate the gradient of the direc- 
tionality W with respect to the expansion coefficients 
{Pm,n}- Due to the large number of involved parame- 
ters, it is crucial to devise a very fast and sufficiently 
accurate scheme to evaluate the gradient, which task can 
be solved with the help of the TCF state basis and its 
self-orthogonality property [47]. 

In Fig. 3d we display the progression of the direction- 
ality V of the first TLM during the optimization proce- 
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FIG. 4. (Color online) Results from the optimization of the 
random laser for emitting into two prescribed directions (65° 
and 190°) simultaneously (the same disorder is used as in 
Fig. 3). Shown are (a) the optimized mode profile |^i(x)| 
inside the random laser, (b) the corresponding pump profile 
resulting from the optimization algorithm and (c) the far-field 
profile featuring a convergence degree of V = 96.12%. The 
inset in (c) shows the near field of the laser mode, just outside 
the pumped disk. Note the very good agreement which we 
find in the main panel of (c) between the optimized emission 
pattern (blue) and the targeted double-peak emission profile 
G(<p) (red). 



dure which leads to a rather quick convergence (in 37 
iterative steps) from an initially random starting point 
to a nearly optimal mode with the desired directional- 
ity. The directionality does not increase monotonically 
as both the gradient of the directionality is only approxi- 
mate and the first TLM can switch between modes. Note 
also that typically several random starting conditions for 
the iteration are necessary due to the existence of local 
maxima of T>({f3 m ^ n }) in the multi-dimensional param- 
eter space of {Pm,n} where the local optimization may 
get stuck. More sophisticated methods to overcome local 
maxima could certainly be applied here. Since, however, 
we already do reach excellent results with V > 96% with 
typically not more than altogether 20 initial guesses, such 
methods are not compulsory. The mode which matches 
the emission pattern of the target function G(cp) best and 
the corresponding pump profile are shown in Fig. 3. This 
encouraging result is directly reflected by the shape of the 
corresponding lasing mode shown in Fig. 3a, c which dis- 
plays the desired narrow, beam-like shape outside of the 
disk boundary. We emphasize, in particular, that the di- 
rectionality can here be much better controlled than in 
strongly disordered media. This is best observed by com- 
paring Figs. 2 and 3: Although the mode which agrees 
best with the desired emission profile G(cp) is selected in 
the strongly scattering medium (Fig. 2), a clear deviation 
from the desired target function G(jp) can be observed. 
In fact, for disorder configurations different than the one 



in Fig. 2 we observed much larger deviations. For weakly 
scattering media, on the contrary, the applied pump pro- 
file is optimized for each specific distribution of scatterers 
such that a directionality of D > 95% is generally guar- 
anteed. 

Our optimization procedure can be pushed even fur- 
ther in the sense that the reference emission pattern G(cp) 
is not limited to a function with a single narrow peak, but 
instead can be chosen almost arbitrarily. As an example 
we have also optimized the pump profile such that the 
random laser emits into two directions simultaneously, 
see Fig. 4. This result shows the versatility of our ap- 
proach and manifests that the emission of a weakly scat- 
tering random laser can be arbitrarily shaped using a 
carefully optimized pump profile. Combining this tech- 
nique together with the frequency selection mechanism 
presented in [31, 32] suggests that random lasers can be 
controlled on an unprecedented scale. 

We emphasize that the pump profiles which give rise 
to these uni- or multidirectional lasing states are nei- 
ther very focused (as in [33, 49]), nor strip-shaped (as in 
[51, 52]), nor in any other way evidently self-explanatory. 
Instead, the obtained profiles are structured on a scale 
comparable to the lasing wavelength - a situation that is 
routinely created in laboratory experiments [14-17, 32]. 
The pump optimization procedure presented above can 
be directly mapped to a corresponding experiment us- 
ing, e.g., a spatial light modulator for shaping the pump 
beam (see Fig. 1). In this case the optimization can be 
organized through a feedback loop between the measured 
emission and the applied pump profile, as implemented 
very recently in a similar setting [32]. Realizing such 
ideas in the experiment would not only be interesting 
for fundamental research on light profile patterning in a 
disordered medium, but would constitute a major step 
forward in bringing random lasers closer to practical ap- 
plications where specific emission profiles are required. 
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